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Abstract 


We study the excited B mesons’ contributions to the coupled-channel effects under the framework 
of the 3P) model for the bottomonium. Contrary to what has been widely accepted, the contributions 
of P wave B mesons are generally the largest, and to some extent, this result is independent of the 
potential parameters. We also push the calculation beyond B(1P) and carefully analyze the contribu- 
tions of B(2S). A form factor is a key ingredient to suppress the contributions of B(2S) for low lying 
[T= bottomonia. However, this suppression mechanism is not efficient for highly excited bottomonia, such 
i as Y(59) and Y(69). We give explanations why this difficulty happens to the 3P) model and suggest 
a analyzing the flux-tube breaking model for the full calculation of coupled-channel effects. 


1 Introduction 


Heavy quarkonium is a multiscale system covering all regimes of quantum chromodynamics (QCD), which 
make it an ideal place to study strong interactions [1]. Despite the success of QCD in the high energy 
== region, due to asymptotic freedom, the nonperturbative effect dominates at low energies and brings 
problems to perturbative calculation. One tool to study this non-perturbative effect is lattice QCD. 
However, due to its huge calculation work, it is still unable to calculate all the physical quantities with 
the current computation power. Another important approach is to develop various phenomenological 
models. Among these phenomenological models, the quark model is a prominent one. Under the quark 
model framework, various types of interactions have been suggested by various groups, and they have 
achieved many impressive successes (see e.g. Refs. [2-6]). However, these potential models cannot be 
the whole story. One important missing ingredient is the mechanism to generate quark-antiquark pairs 
which enlarge the Fock space of the initial state, i.e., the initial state contains multiquark components. 
These multiquark components will change the Hamiltonian of the potential model, causing mass shift and 
mixing between states with the same quantum numbers or directly contributing to an open channel strong 
decay if the initial state is above the corresponding threshold. These consequences can be summarized as 
unquenched effects or coupled-channel effects. 

Through various approaches, such as the °P) model [7-11], flux-tube breaking model [12-16], micro- 
scopic decay models [4, 17-19], S matrix analysis [20-23], coupled-channel effects are extensively studied 
in many literatures (e.g., Refs. [24—29]) 
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Despite these pioneering works, we find that at least two factors have the potential to jeopardize the 
calculations of the coupled-channel effects. One is the widely used simple harmonic oscillator (SHO) wave 
function which approximates the realistic wave function in the wave function overlap integration. As 
already pointed out in our previous work [30], the SHO approximation is not good, especially for states 
near thresholds, and it is essential to treat the wave functions precisely. 

Another is the assumption that the contribution of the excited meson loops is negligible. Take the 
bottomonium, which is the system we study in this paper, as an example. Two reasons may explain 
why this approximation is widely used. One is that the threshold corresponding to the excited B mesons 
is higher than the ground state BB threshold; thus the coupled-channel effects is expected to be small. 
Another is that the calculations of the excited B mesons’ contributions are more complicated, an effective 
method to do the calculation is still not widely known. 

As explained in [30], the Gaussian expansion method plus the techniques of transformation between 
the Cartesian and spherical basis can, in principle, do the sophisticated calculations. Another observation 
is that the first excited P wave B meson is only around 450MeV heavier than the ground state B meson. 
This mass difference is less than one tenth of the B mesons’ mass. Given the fact that the quantum 
numbers of Bı are different from the ground state B mesons, and the coupled-channel effects do rely on 
the quantum numbers, the suppression purely originates from the larger mass may be not as large as what 
has been taken for granted. 

The sum of all the intermediate meson loops is more than just a calculation challenge, it may also 
lead to profound physics. In the light sector, a series of works by Geiger et al. [31-33] show us that even 
though different intermediate meson loops contribute to the breaking of Okubo-Zweig-lizuka (OZI) rule, 
under some simplifications (such as neglecting the mass difference in the denominator), they contribute 
destructively, leaving us a perfect OZI rule. If one only sums over some of the meson loops, one may leave 
with a wrong conclusion. 

One should not confuse their calculation with what we are going to do in this work. They studied the 
flavor-changing process, such as uii + virtual meson pairs > dd. However, we mainly focus on the mass 
shift and their cancellation does not happen to our case. For states below the threshold, the intermediate 
loops always contribute a negative mass shift. Even though for the above threshold case, where the mass 
shift may add destructively, the mass difference does matter in the real calculation. 

In Ref. [15], under the flux-tube breaking model, Geiger and Isgur also showed that if one adds up 
all the intermediate states, the mass shift caused by coupled-channel loops does not converge when two 
assumptions are adopted; i) the meson wave function is the SHO wave function, and ii) the string length 
between the generated quark pairs is zero. In the ĉP) model, this conclusion needs to be checked or 
recalculated, because the SHO approximation is far from true, which was already stressed before and the 
vertex of the quark generation is different from the flux-tube breaking model. 

There are also some studies of the excited meson loops for heavy quarkonium. In the charmonium 
sector, a lattice calculation in Ref. [34] shows that 7. and J/w has small but non-negligible components 
of D,D*. In bottomonium sector, loops involving a Bı meson are proven to be critical to explain the 
large breaking of the heavy quark spin symmetry of Y (10860) [35]. However, these studies focus on some 
specific states, and systematic studies of the excited meson loops is still missing. 

To summarize, it remains to be answered or clarified whether the ground state approximation is good 
or not, and the general properties of the excited meson loops still need to be systematically exploited. In 
this work, we try to answer these questions under the framework of the 3P) model. 

This paper is organized as follows. In Sec. 2, we briefly describe the ingredients of the coupled-channel 
effects and the calculation techniques, including the 3P) model, the Cornell potential model, the Gaussian 
expansion method, and the transformations between the spherical and Cartesian basis. Sec. 3 is devoted 
to the results of the higher excited B mesons up to B(2S), and we explain the necessity of the form factor 
and the limitation of the 3P) model. Finally, we give a short summary in Sec. 4. 
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2 Theoretical Framework 


The calculation methods and tools are described in details in Ref. [30]; thus, we only sketch the key steps 
in this section. 


Figure 1: Sketch of coupled-channel effects in the °P) model for the bottomonium. i and f, respectively, 
denote the initial and final states with the same JPC and BB stands for all possible B meson pairs, 
including the excited B mesons. 


The coupled-channel effects in the 3P) model [7-9] can be described by Fig. 1. In this model, the 
generated quark pairs have a vacuum quantum number JPC = 0+*. In the notation of STI Ly, it is 3P) 
that explains the model’s name. 

The Hamiltonian to generate the quark-antiquark pairs is expressed as 


A; = amay | Body (1) 


where mz, is the produced quark mass andy is the dimensionless coupling constant. Since the probability to 
generate heavier quarks is suppressed, we use the effective strength 7, = at in the following calculation, 
where mg = My, = Ma is the constituent quark mass of the up (or down) quark and ms is the strange 
quark mass. 

As sketched in Fig. 1, the experimentally observed state should be a mixture of a pure quarkonium 
state |o) and a B meson continuum state |BC;p). The experimentally observed state |A) should be 
expressed as 


A) = alo) +> | dpeno()|BC:p), (2) 
BC 


which is the eigenstate of the full Hamiltonian H defined in Eq. 6. co and cgc stand for the normalization 
constants of the bare state and B meson continuum, respectively. 

In order to work out |A), as a first step, one has to solve the wave function |) for the heavy 
quarkonium. In this work, it is obtained by solving the Schrödinger equation 


nies (2m $ z prea V(r) Fo) = alo} (3) 


where mp, and Mo represent the mass of b quark and the bare mass of the bottomonium, respectively. In 
the above equation, V(r) is the well-known Cornell potential [4, 17] 


V(r) = 22 4 arte, (4) 
3r 


where a,A and c stand for the strength of color Coulomb potential, the strength of linear confinement 
and mass renormalization, respectively. V,(r) stands for the spin dependent interactions which restores 
the hyperfine or fine structures of the bottomonium, 


2a À 327 da (o (Sp F)(SE-7) 
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where L denotes relative orbital angular momentum, S= S, + S; i is the total spin of the b quark pairs and 
mp is the b quark mass. ô(r) is the smeared delta function and is written as ô(r) = (o/,/m)3e-?'" [19,36]. 
We treat the spin dependent term as a perturbation, and the spatial wave functions are obtained by 
solving Schrödinger equation numerically using Numerov’s method [37]. 

Combine the Cornell potential and the dynamics of quark pair generation, we get the full Hamiltonian, 


H = Ho + Ego + Hr, (6) 


where Epc = y m? +p? + v me, + p°. The H; term in Eq. 6 is mainly responsible for the mass shift. As 
the name tells us, it is naturally defined as 


AM = M — Mo (7) 
and it can be obtained by solving the integral equation 
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(8) 


Note that the ze term is added to handle the situation when m4 > mpg + moc. In this case, AM will pick 
up an imaginary part 


D a (BC; Pa|Hil|¥o) |’, (9) 


which is equal to one half of the the decay width. Pg and Eg denote the momentum and energy of B 
meson, respectively. The wave function overlap integration lies in the term 


(BC; Ps|Hrlyo) = > | Pkdo(k + Pa)op(k + ePp)do(k + xPp)|k|Y{" (Gz, oz) (10) 
polarization 


where £z = mg/ (Ma + mq) and Mm, and mg denote the b quark and the light quark mass, respectively. 
Once M is solved, the coefficient of different components can be worked out either. For states below 
the threshold, the probability of the bb component is expressed as 


Ls E 
Pg= a= (14 >. dp er) (11) 


BCLS 
where |M°5|? is represented as 
MESI = | ds (BC; Pal tile)? (12) 
The main calculation work lies in the Eq. 10. In order to evaluate it precisely, we use the Gaussian 


expansion method (GEM) [38] and the transformation between the spherical and Cartesian basis to make 
the GEM automatic applicable to the excited B meson. 


3 Results and Discussions 


3.1 B,’s Contributions with Traditional °P) Model 


In this work, we do not intend to reproduce the spectrum and decay widths of the bottomonium family, 
but to study the general contributions of the excited B meson to coupled-channel effects. 

There are two simplifications in our calculation. Firstly, since the constituent quark mass of u, d quarks 
are set to be the same to get the wave functions, we further make the approximation m(B°) ~ m(B* 
Notwithstanding this simplification, there are still 42 channels to be calculated, i.e., B im B 7 ,B B,(1P) 
and Bis) (1P) Bis) (1P). In principle, one can still treat their mass precisely, however, based on our expe- 
riences, this simplification is quite precise, and it saves a lot of the calculation work. 

Second simplification comes from Bı meson multiplets. From the perspective of the quark model, 
there are four 1P wave B mesons. Nevertheless, the predicted Bo(1P) and one of the B,(1P) are still not 
experimentally observed. As far as the coupled-channel is concerned, the missing of these two states is 
not a real problem. Because the mass of b quark is very large, it is safe to assume that the heavy quark 
spin symmetry does not break. Under this limit, it is reasonable to set the mass of Bo(*?Po) and missing 
partner of B,(5721) to be the same as B,(5721). There is a mixing between B(?P,) and B('P;) to form 
the experimentally observed B,(5721) [39]. However, since we sum up all the possible combinations of B 
mesons, mixing between the two 1*B mesons can also be neglected. As far the mass of the intermediate 
states are well measured, we use the value from the Particle Data Group [39]. 

The coupled-channels are much more involved than the case, which only includes the ground state B 
mesons. For illustration purposes, we sum up the contributions of the states in the same multiplet and 
u, d,s flavors and their charge conjugate partners, classifying the channels into three groups (see Fig 2), 
e.g., B(1S) stands for ground state B, B*, B,, and BS mesons. Also, in the rest of this paper, we define 
AM(i,j) to represent the sum of the mass shift due to all the possible combinations of i and j wave B 
mesons. 

The parameters are given in Table 1, which are the same as our previous work [30]. 


SS 


a = 0.34 A = 0.22GeV? c = 0.435GeV 
mp = 4.5GeV my, = Ma = 0.33GeV m, = 0.5GeV 
o = 3.838GeV y = 0.205 


Table 1: The parameters used in our calculation. 


The mass shift AM of both B(1S) and B(1P) mesons are depicted in Fig. 2. In the rest of this 
section, we are going to discuss some interesting structures of the plots. 

When the states are far below the threshold, the mass shift is always negative [as indicated in Eq. 8], 
and closer to the threshold means an increase of |AM|. When the mass goes higher, the magnitude of 
AM is no longer monotonically increasing but oscillating, which reflects the node structure of the wave 
functions. This phenomenon happens to both cases whether Bı mesons’ contributions are considered or 
not. 

A specific example would help to clarify this point. For Y(10860), which we treat it as T(5S), 
AM(15S,1P) is generally the largest but there is a big dip around 11.14 GeV with parameters in Tab. 1; 
however, in Ref. [40]’s parameters, this channel’s contribution always grows. In these cases, it is difficult 
to make a solid conclusion about the spectrum behavior, i.e., the spectrum is sensitive to the potential 
parameters. This conclusion also agrees with Ref. [41]. 

The most notable impact of Bı family is the unexpected large contributions to the mass shift, and 
generally, AM(1S,1P) is the largest. Even for Y(1S), which is far below the threshold, AM(1S,1P) + 
AM (1P,1P) is 14 times larger than AM(15, 1S). 

Compared with the case that only considers the ground states B mesons contributions, threshold effect 


is more clearly reflected when the excited B mesons is included, because open channels due to different 
multiplets are well separated. 
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Figure 2: —AM of different coupled-channels for Y(nS). The results with parameters in Tab. 1 and the 
results recalculated with parameters of Ref. [40] are represented by thick and thin curves, respectively. 
AM(1S,1S), AM(1S,1P), and AM(1P,1P) are represented by black dot-dashed, red dashed and blue 
solid curves, respectively. 
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Figure 3: —AM of different coupled-channels for some selected representatives of the tSo, 1P}, 3 P3, 3 D1 
families. The results with parameters in Tab. 1, and results recalculated with parameters of Ref. [40] 
are represented by thick and thin curves, respectively. AM(15,15), AM(1S,1P), and AM(1P,1P) are 
represented by black dot-dashed, red dashed, and blue solid curves, respectively. Here we omit the results 
of (nS), hy(nP), and x,;(nP) when n < 2, because these results are quite similar to those of T(1S) or 
T(28). 


The first open bottom threshold is 2mg œ~ 10.56GeV. When the mass approaches to this value, 
the denominator M — Ego approaches to 0 in Eq. 8; thus, there is an enhancement of AM from the 
B(1S)B(1S) channel. From the (39) case in Fig. 2, one can clearly see the sharp increase of AM(15, 1S) 
between 10.55 and 10.6GeV, while the slop of the other channels does not change much because their 
threshold is 450MeV larger. 

The threshold effect is broken in some degree by the nodes of the wave functions. The peaks and 
valleys of the wave function are more likely to cancel with each other for higher excited states, leaving a 
relatively small slope of AM. One can compare the result of T (4S) with Y(3S) verify this conclusion. 

We need to stress that 3S1 is not the only family who couples strongly with BB(1P) loop. As shown 
in Fig 3, all the families of 1 Sọ, P1, 3 P3,3 Dı share these general properties. 

An direct consequence of the large AM(15S,1P) is that the parameters only considering ground state 
B mesons to reproduce the experimental data are somewhat incomplete or even misleading. 

This conclusion is independent of the parameters to some degree, since the results are based on two 
different sets of parameters (Refs. [30] and [40]), both give large AM(1S,1P). 

It seems feasible that to fit the spectrum, one may tune the parameters, such as the mass renormal- 
ization term c in Eq. 4, to absorb this unexpected large mass shift, or take a further step, weaken the 
confine potential to reflect the fact that closer to the threshold generally means more suppression of the 
mass. 

This could be a possible way out, if someone only focuses on the spectrum. However, as pointed out 
in our previous work [30], the mass shift only reveals one aspect of the coupled-channel effects. 

To be specific, the big contribution of B(1P) not only brings us the large mass shift AM, but also 
the large fraction of meson pairs (or equivalently, small P,z;). This renormalization of the wave function 
cannot be accommodated in the framework of the potential model. 


3.2 Beyond Bı 


The large contribution of B,; mesons reminds us that it is essential to do the ab initio calculations of 
coupled-channel effects. However, before really carrying out these calculations, one may naturally ask this 
question, “how large are the contributions of the states beyond Bı or how to evaluate all the intermediate 
meson loops?” 

It is very difficult to offer a complete answer to this question. Nevertheless, it turns out that we can 
still estimate the contributions up to B(2S). Furthermore, by analyzing the mass shift behavior carefully, 
we can draw some model independent conclusions. 

First of all, we are trying to precisely evaluate the loops. Of course, since the B(2S) mass and its 
wave functions are not known, the calculation has to be model dependent, which is inevitable. 

So far, the B(2S) meson has not been experimentally well determined and its mass is model dependent. 
For our purposes, it is sufficient to take an average of several theoretical estimations in Refs. [42—46], and 
for the consistency treatment of the wave functions, the B(2S) wave functions are still deduced from our 
parameters. 


We define the ratio 
AM (2S, 2S) + AM(2S,1P) + AM(2S, 1S) 


AM(15, 15) + AM(15, LP) + AM(1P, 1P) 


to tell how large are the contributions of B(2S). 

As will be explained below, in the calculation when B(2S) are involved, we adopt the Gaussian form 
factor which reveals the size of the generated quark pairs (see, e.g., Ref. [47-49]). This form factor modifies 
(BC; Pg|Hr|Wo) to be 


R:= 


(13) 


(BC; Pa\Hilo) = XO | Pk” polk + Pao (k + cPp)dh(k + 2Pp)|KIYi"(0z,2), (14) 


polarization 


In this work, we fit r to be 0.408 fm, which minimizes R of Y(1S). This value slightly larger than the 
value r = 0.335 fm used in Ref. [49]. 

From Fig. 4, one can clearly observe that the form factor is crucial to suppress the B(2')’s contribu- 
tions. Without it, B(2S) can contribute an additional 25% of the AM. This deviation is not negligible 
if one wants to make a precise fit of the spectrum. However, with a form factor r = 0.408 fm, the AM 
of B(2S) is suppressed to less than 5%. The previous question is partly answered, the form factor is 
indispensable if one needs to suppress the contributions of B(2S). However, with the increase of n in 
Y(nS) family, the form factor works less and less efficiently to suppress R. 
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Figure 4: R [defined in Eq. 13] of the Y family. Results corresponding to r = 0 (no form factors), 
r = 0.335fm (taken from [49]), and r = 0.408fm (our best fit) are denoted by black dot-dashed, blue 
dashed, and red solid curves, respectively. 


Another effect of the form factor is that it adds more peaks to the oscillation of AM. This result 
has a pictorial explanation. Since the form factor effectively serves as a cutoff, this term will generally 
suppress every channel’s contribution. The key point is that different channels are suppressed by a different 
magnitude. With the increase of the nodes of the wave functions of bottomonia and B mesons, at some 
specific energy point, B(2S)’s relative contributions may increase. 
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Figure 5: -AM of Y(5S) and Y(6S) with different form factors. Solid curves represent the denominator 
of Eq. 13, i.e., AM(1S,1S) + M(1S,1P) + M(1P,1P), and dashed curves represent the numerator, i.e. 
AM (28S, 2S) + AM(2S,1P) + AM(2S, 1S). AM corresponding to r = 0, r = 0.335fm, and r = 0.408fm 
are denoted by black, blue, and red curves, respectively. The value of the -AM cannot compare within 
different form factors because the 3 P) coupling constants are not fit to reproduce experimental data. 


Clearly the sharp peaks of T(5S) and Y(6S) deserves a special analysis. Both of them locate around 


10 


11.18GeV. Given the fact that the wave functions of Y(5S) and Y(6S) are quite different, yet they both 
have the same sharp peaks at the same energy points, it turns out to be more than just a coincidence. 

R is convenient to estimate B(2S)’s relative contributions; however, it cannot tell whether the peak 
is due to the suppression of the denominator or the increase of the numerator. The complete information 
is encoded in AM itself. In Fig. 5, we show AMs for Y(5S) and Y(6S) with different form factors. 

From Fig. 5, one can see that the enhancement of R around 11.18 GeV is in fact the bigger suppression 
of AM(1S,1S) + AM(1S,1P) + AM(1P,1P), [or more precisely, the suppression of AM(1S,1P)]. We 
have scanned r in quite a wide range 0.17 — 1.08 fm, which all show that the form factor does not shift 
the positions of peaks or valleys of AM. Around this region, the form factor fails to suppress the B(2S)’s 
contribution is somewhat a coincidence, because this suppression does not happen to T(5S) and Y(6S) 
with Ref. [40]’s parameters; their AM(1S,1P) always goes up. 

Even in the low energy range where all the coupled-channels contribute a negative AM, no matter 
whether the form factor is considered or not, the B(2S) still contribute around 20% of the total mass 
shift for (5S) and Y(6S). In other words, the form factor fails to suppress B(2S)’s contributions. This 
strongly indicates that only adding the Gaussian form factor is not adequate to result in a complete 
calculation of the coupled-channel effects under the 3P) framework. 

One may argue that the coupled channel effect can be absorbed into the smeared potential which 
returns wave functions broader than those of the Cornell model. These broad wave functions indeed 
somewhat suppress the higher excited states contributions. However, based on our experience, changing 
the size of the wave functions is also not able to solve the convergence issue. 

We conclude this section with some comments on the °P) model. In the classical Py model, the 
suppression due to the natural cutoff from wave functions and the increase of the denominator of Eq. 8 
are too weak. Even with the remedy of the Gaussian form factor, if one does not choose the cutoff r 
carefully, one will enhance instead of suppressing the B(2S) contribution. Compared with other models, 
such as the microscopic decay model or the flux-tube breaking model, which have rich microscopic details 
of the decay vertices, the 3° P) model replaces these fine structures as an overall coupling constant y. This 
approximation may be not appropriate (see e.g. Ref. [50]). In our case, we show that this approximation 
leads to the bad convergence of the sum of the excited meson loops for highly excited bottomonia. 

The rich structure of the quark pair generation vertices may help to solve the convergence issue. For 
example, the Hamiltonian of the flux-tube breaking model suppresses the generation of the farther quarks, 
one will expect that it is harder to generate excited intermediate states. 


Figure 6: Sketch of coupled-channel effects in the flux-tube breaking model. Red lines stand for the 
flux-tube. The regions outside the flux-tube are shaved off, i.e., excited states’ contributions in the loops 
are suppressed. 


As illustrated in Fig. 6, contributions from the regions that are farther from the flux-tube line can be 
dropped, as a result, the excited mesons’ contribution is naturally suppressed. Without this dynamical 
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suppression, it is more difficult to suppress B(2S)’s contributions even with the modified version of the 
3 Py model, where the form factor is added. 


4 Summary and Outlook 


In this paper, under the 3° P) model framework, we explicitly calculated the excited B mesons’ contributions 
to the coupled-channel effects for the bottomonium. We reveal the fact that compared to the ground 
state B mesons, contributions from B(1P) mesons are generally the largest. Up to this partial wave, it is 
necessary to do the ab initio calculations of the coupled-channel effects. 

When we push the calculation beyond B(1P), we find some fundamental difficulties of the 3P) model. 
Even with the carefully chosen form factor, it still cannot efficiently suppress intermediate state contri- 
butions of higher partial waves. Since we do not fit our parameters with experimental data, and we have 
exploited several different sets of parameters, we have enough reasons to believe that the difficulties are 
independent of the wave functions or the potential models. 

We suggest that an efficient suppression mechanism such as a dynamically suppression is needed to 
evaluate the coupled-channel effects. How to effectively sum up all the intermediate loops of coupled- 
channels still remains to be an open issue. 
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